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Abstract. We present UBVRIz' photometry of the gra- 
vitational lens candidate CLASS B1152+119 obtained 
with the Nordic Optical Telescope. The two QSO com- 
ponents are resolved in the B, V, R, I and z' bands con- 
firming the lensing nature of the system. The z = 0.44 lens 
galaxy is clearly detected in B, R, I and z' and its posi- 
tion is found to be almost coincident with the faint QSO 
image which is heavily extincted (relative to the brighter 
QSO image) by dust in the lens galaxy. The extinction 
curve of the lens galaxy derived from the relative pho- 
tometry is well fitted by a Galactic extinction law with 
1-3 <Rv< 2.0 and E{B - V) « 1. From a simple model 
of the system we predict a time delay of ~ 60 days. 

Key words: quasars: individual: B1152+199 - galaxies: 
ISM - ISM: dust, extinction - ISM: evolution - gravita- 
tional lensing 
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Dust is ubiquitous in the Universe and is responsible for 
extinction of the light from distant sources. The type of 
dust (composition, grain sizes and shape) determines the 
amount of extinction as a function of wavelength - the 
extinction curve A\. Extinction curves almost certainly 
evolve with redshift since the metallicity, elemental abun- 
dance ratios, mean star-formation rate, and energy-injec- 
tion rates that determine the structure and evolution of 
the dust are all strong functions of redshift. The extinction 
in our Galaxy has been well studied (Cardelli et al. 1989), 
but very little is known about extinction in other galax- 
ies, especially at higher redshifts. This is unfortunate as 
dust plays an increasingly important role in contemporary 
cosmology. 

Recently the use of multiply imaged QSOs has been 
explored as a means of inferring differential extinction 
curves for distant galaxies (Nadeau et al. 1991; Malho- 
tra et al. |1997|; Falco et al. [19990. The idea behind this 



* Nordic Optical Telescope is operated on the island of La 
Palma jointly by Denmark, Finland, Iceland, Norway, and Swe- 
den, in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. 



approach is that a lensed QSO can be used as a 'stan- 
dard beacon', shining through different paths of the lens- 
ing galaxy. If one of the images suffers negligible extinc- 
tion, or if the extinction curve is the same throughout the 
galaxy, then the relative intensity ratios of any two im- 
ages as a function of wavelength is a direct measure of the 
extinction curve. 

As a part of the CLASS radio survey, Myers et al. 
( 1999| ) reported the discovery of a new gravitational lens 
candidate, CLASS BT152+199, with an image separation 
of 9 = l'/56 and a flux ratio of 3.03 ± 0.03 at 8.46 GHz. 
Spectra obtained with the Keck II telescope revealed a 
background quasar at z = 1.019 and a foreground galaxy 
at z = 0.439. Optical follow-up observations revealed a 
bright unresolved object (g ~ 16.5 and i ~ 16.6) located 
at the radio position, presumably the brighter of the two 
lensed images. The non-detection of the weak image in 
g and i is indicative of a large color difference between 
the two images. If B1152+199 is indeed a gravitationally 
lensed system this could be attributed to differential ex- 
tinction caused by dust in the lens, making BT152+199 a 
well-suited system for studying the extinction law of the 
lensing galaxy. 



2. Observations and reduction 

We have performed multiwavelength U BVRIz' photome- 
try on data of B1152+199 obtained with the ALFOSC in- 
strument at the 2.56-m Nordic Optical Telescope (NOT). 
The ALFOSC instrument is equipped with a Loral/Lesser 
CCD with an approximate gain of q = le~/ADU and 
RON = 6 e" /pixel and a pixel scale of O'.'188/pixel. The 
observations were performed within the span of a few days 
to avoid time delay and variability effects. The observing 
log is given in Table |l|. The conditions were photometric 
with average seeing ranging from (f.'8 in the z' band to l'/3 
in the U band. The CCD frames were reduced (overscan 
subtracted, trimmed and flatfielded) using an IRAF-based 
pipeline. [J] 
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Table 1. Observing log giving the dates of observation 
and the total exposure times for each night 



B1152+199 



On June 25 1999 UT, we also obtained photometric 
UBVRI observations of standard stars in M92. Using 
these and calibrated photometry by Davis (private com- 
munication) we fitted the transformation equations to 19 
stars in M92 to derive the color, extinction and zero-point 
terms for the calibration to the Johnson-Kron-Cousins 
system. 



PSF star 



3. Deconvolution of B1152+199 

The two QSO images together with the lens galaxy were 
deconvolved with the MCS deconvolution algorithm (Ma- 
gain et al. 1998] ). This method is based on the principle 
that the resolution of a deconvolved image must be com- 
patible with its sampling, which is limited by the Nyquist 
frequency. In order to improve the sampling of the images 
we adopted a pixel size in the deconvolved images as half 
the pixel size of the o riginal frames (C/188/2 = 0'.'094). 
(See Burud et al. 1998| for more details on the implemen- 
tation of the MCS algorithm). 

In order to successfully deconvolve the images, the 
point-spread function (PSF) must be bright and located 
close to B1152+199 to minimize the photon noise and 
to avoid PSF variation as a function of the position on 
the chip. We used the bright star 13" South West of 
B1152+199 as a PSF star (Fig. [I]). Careful examination 
revealed that this star had a faint red companion. By de- 
convolving this star with another PSF in the field we were 
able to resolve the red companion star and subtract it, be- 
fore using the bright star as a model for the PSF. We com- 
pared the star to other stars in the field to check that the 
companion star was correctly subtracted. The best fitting 
positions of the weak QSO image relative to the strong 
QSO image returned from the MCS algorithm is given 
in Table @. These yield a smaller separation that expected 



from the radio positions (Myers et al. 1999). The deconvo- 
lution code is very sensitive to PSF variations, and small 
residuals from the subtraction of the red companion could 
explain the discrepancies between the fitted positions and 
the radio position of the weak QSO image. 

In order to improve the photometry we fixed the po- 
sition of the weak QSO image relative to the bright QSO 
image to the radio coordinates. The quality of the re- 
sults was checked from the residual maps, as explained 
by Courbin et al. ( [1998 ). During the deconvolution the 
lens galaxy can be modeled numerically or by an an- 
alytical function (e.g., de Vaucouleurs, exponential-disk 



Fig. 1. 160" x 140" finding chart for B1152+119. North is 
up, East is to the left. The star used to model the point- 
spread function is indicated 

law). The lens galaxy in B1152+199 was well fitted by the 
numerical model and an exponential disk law, whereas 
a de Vaucouleurs model gave significantly bad residu- 
als. The best fitting exponential disk law was centered 
at (Aa, AS) = (0"29 ± 0"02, -0'.'89 ± 0"04) relative to the 
bright QSO image. 

Both components of B1152+199 were resolved in the 
V, R, I and z' bands (Fig. |). Following Myers et al. fll999| ) 
we will refer to the bright component as "A" and the weak 
component as "B" . In the U and B bands we failed to de- 
tect the B image due to the increasing size of the seeing 
disc with shorter wavelength and the large differential ex- 
tinction of the system (see Sect. ||). The lens galaxy was 
detected in the R. I and z' bands. We performed aperture 
photometry of the galaxy on the deconvolved images with 
an aperture diameter of 2"26. The U, B, V, R and / band 
magnitudes of the A and B images and the galaxy (G) 
were calibrated (Table [5j) whereas no standard photome- 
try was available in the z' band. 

Additional extended signal was observed just South 
East of the B image (notably in the I band). We discuss 
this further in Sect. |]. 



4. PSF photometry of B1152+199 

In addition to the deconvolution we performed PSF pho- 
tometry on the system with the IRAF DAOPHOT pack- 



age (Stetson 1987). As PSF star we used the same star 



as in the deconvolution, with the faint red companion re- 
moved with the IRAF task substar. Again we fixed the 
position of the B image relative to the A image on the the 
radio position to ensure the correct separation. The results 
from this analysis is not as robust as the results from the 
deconvolution since it does not simultaneously fit a model 
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Fig. 2. Deconvolved images of B1152+199 in the (from 
left to right) z', R and / bands. The resolution (FWHM) 
is 0'.'28 in the R band image and 0'.'19 in the / and z' 
band images. The bright QSO component is labeled "A" 
and the faint component is labeled "B" 




for the galaxy, but the differences between the results of 
the two approaches is indicative of the magnitude of the 
systematic photometric uncertainties. 

The two QSO images were resolved in the B, V, R, 
I and z' bands. The results from the PSF photometry is 
given in Table |[ After subtracting the scaled PSFs we 
detected the lensing galaxy in the B, R, I and z 1 band 
residual images (Fig. ||). The / and z' band images indi- 
cates that the extended flux South East of the B image, 
detected in the deconvolved images, stems from the lensing 
galaxy. The light from the B component passes through 
the galaxy close to its center, and in the deconvolution 
process some of the light from the galaxy may have been 
attributed to the B image. 

To test the results from the deconvolution we per- 
formed aperture photometry with the SExtracter soft- 



ware (Bertin & Arnouts 1996) on the images cleaned with 
DAOPHOT. The residual images in Fig. || show artifacts 
of the PSF subtraction (seen as circular dark spots). These 
affect the aperture photometry which should therefore be 
considered rather uncertain. In the / and B bands we 
obtain / = 19.43 ± 0.03 and B = 20.24 ± 0.02 with an 
aperture diameter of 2'.' 26. In the R band the signal is 
rather weak, and thus in the outer parts dominated by 
the PSF subtraction artifacts. We can however obtain a 
color by comparing the R and / flux within an aperture 
of T/13. This yields R-I = 0.80 ± 0.05. The color of the 
galaxy from the deconvolution (R-I)mcs — 0.43±0.06 or 
from PSF photometry (R — I)psf = 0.80 ±0.05 are consis- 
tent with the lensing galaxy being a irregular or late type 
galaxy at z = 0.439 (Fukugita et al. 1995). By comparing 
the photometric results from the two approaches we esti- 
mate the magnitude of the systematic photometric errors 
due to the limited quality of the data to be in the order of 
0.4 mag for galaxy magnitudes and 0.2 mag for the QSO 
magnitudes. 




Fig. 3. PSF-subtr acted z' (top left), / (top right), R (bot- 
tom left), and B (bottom right) images 
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Table 2. Columns 2-4 give the flux ratios and astrome- 
try for B1152+199 from the deconvolution. Aa and A<5 
are the fitted positions of the B image relative to the 
A image. The rms uncertainties in the positions are ap- 
proximately 0'.'05 in Aa and O'.'Ol in AS. The positions 
agree reasonably well with the radio positions (Aa, AS) = 
(0'.'935 ± 0'.'005, -17248 ± 0'.'005), although the separation 
is slightly smaller. This is a well-known effect due to the 
entanglement of the B image and galaxy flux. Column 5 
gives the flux ratios from the PSF photometry 
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Table 3. Photometry (mag) of image A, B and the lens- 
ing galaxy. The uncertainties are la standard deviations, 
including standard-star zero point, calibration and pho- 
tometry uncertainties 
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5. Galaxy extinction 

Significant differential extinction of the B image relative to 
the A image was observed (Tables | and |). The A image, 
however, does not appear to be significantly reddened. In 
what follows we assume that the A image is unaffected 
by extinction by dust in the lensing galaxy. The Galactic 
extinction curve is fairly well known, with A\ well param- 
eterized as a function of a single parameter, Ry, defined 
as the ratio between the V extin ction , Ay, and th e colo r 
index E(B - V) (Cardelli et al. |l989t Fitzpatrick |l999f) . 
This parameterization is valid over the observed Ry val- 
ues ranging from 2.5 to 6, with a mean (canonical) Galac- 
tic value of Ry =3.1. Large Ry give fiat UV extinction 
curves due to large grains whereas small Ry and small 
grains lead to large UV extinctions. 

The variation of the flux ratio with wavelength directly 
measures the extinction curve of the galaxy. We have as- 
sumed that the intrinsic flux ratio is given by the radio 
flux ratio (B-A=1.21). 

The extinction curve of B1152+199 was well fitted by 
a Galactic extinction law, using either the flux ratios from 
the deconvolution or from the PSF photometry. However 
the best fitting values of the parameters Ry and E(B — V) 
differ somewhat. We first discuss the fit to the deconvolu- 
tion data. In the deconvolved U and B band images the 
galaxy and the B image were too weak to be separated. We 
can however obtain upper limits on the B image flux by us- 
ing the total flux of the galaxy and the B image. These are 
B 

gaiaxy+B — 20.37 and C/gaiaxy+B — 20.45. The best fitting 
restframe values of the parameters to the deconvolution 
flux ratios are Ry = 2.04±0.04 and E(B-V) = 0.89±0.02 
(Fig. |J, solid curve). The predicted magnitude of the B im- 
age from the observed A image magnitude and Ry = 2.04 
are Bb = 21.69 and Ub — 23.27. Furthermore the inferred 
-Bgaiaxy+B agree well with the B band magnitude of the 
galaxy inferred from the PSF photometry (Sect. |J). It is 
therefore likely that the upper limits presented above are 
dominated by the galaxy flux rather than flux from the B 
image. 

The PSF photometry detected the B component in the 
B band which gives us another point on the extinction 
curve. The best fitting restframe values of the parameters 
to the PSF photometry flux ratios are Ry — 1.32 ± 0.02 
and E(B - V) = 1.13 ± 0.01 (Fig. |, dashed curve). The 
extinction curves inferred from the two approaches are 
comparable. As discussed in Sect. [| the systematic un- 
certainties in the photometry is of the order 0.2 mag- 
nitudes. This uncertainty is illustrated by the different 
slopes of the extinction curves derived from the two ap- 
proaches. We conclude that the extinction curve of the 
lensing galaxy of Bl 152+199 is fitted by the Galactic ex- 
tinction law and put the following constrains on the pa- 
rameters: 1.3 <R V < 2.1 and 0.9 < E(B -V)< 1.1. 

If the reddening is due to dust extinction in the 
lens galaxy we can derive the column density of neu- 
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Fig. 4. Extinction curve of the lensing galaxy of 
B1152+199, fitted by a Galactic extinction law. The best 
fitting parameters to the deconvolution data (filled sym- 
bols) are Ry = 2.04 ± 0.04 and E(B - V) = 0.89 ± 0.02 
(solid curve). The U and B band points are lower limits 
(see text). The best fitting parameters to the PSF pho- 
tometry data (open symbols) are Ry — 1.32 ± 0.02 and 
E(B — V) = 1.13 ± 0.01 (dashed curve). For comparison 
we have overplotted the Galactic (Ry = 3.1) extinction 
curve (dotted curve) 
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Table 4. Simple lens model parameters for B1152+199 
(cf. Fig. ||). Act, AS are the positions of the lensing galaxy 
relative to the A component. The position of the galaxy 
and the QSO images along with their flux ratios were held 
constant in the fitting procedure, while the time delay was 
calculated from the best fit to the remaining parameters, 
fi = 0.2, A = 0, Hq = 65 km s _1 Mpc -1 were assumed 



tral hydrogen from the observed E{B — V) by assum- 
ing a gas to dust ratio for the galaxy. For the Milky 
Way (N(HI)/E(B- V))„ w = 4.93 x 10 21 cm" 2 mag" 1 , 
(Diplas & Savage 1994 ) while the value for Damped 



Lymana Ab sorbe rs (DLAs) may be up to 20 times larger 
(Fall & Pei |1993|) . With the inferred E(B - V) w 1 for 
the lens galaxy: N(H I)bii52+199 = k x 10 21 cm -2 , where 
5 ^ k < 100 is the parameterization of the range in the 
observed (N(HI)/E(B — V)) in different systems. Thus, 
the lensing galaxy of B1152+199 would qualify as a DLA 
(N(HI)> 2 x 10 20 cm -2 ) regardless of the value of k. 



6. A simple model of B1152+199 

A simple i 
S(i?) = So 



A simple model with a dark halo, represented by 

^"core^cut ( ' ' 



7"cut 



(r 2 orc + J R 2 ) (r 2 nt +i? 2 ), 
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Fig. 5. Critical (solid) and caustic (dashed) curves from 
the model, overlaid the deconvolved R image 

as the dominating mass was fitted to the data (Kneib et 
al. 1996| ). Although the observational constraints available 
from our data were too poor to perform detailed modeling, 
a set of best-fitting parameters could be derived for the 
system (see Table ||). We used the radio positions and flux 
ratios for the two QSO components while we experimented 
with different positions for the lens galaxy. We found that 
the two-image configuration could not be reproduced with 
the galaxy in the position found from the deconvolved im- 
ages. An improved fit was obtained by placing the galaxy 
closer to the B component. The position was chosen to 
be consistent with what is found in the PSF-subtracted 
images. This position means that the B component passes 
through the galaxy while the A component passes outside 
the galaxy, which is consistent with the large observed dif- 
ferential extinction. The model predicts a time delay of 59 
days. In Fig. || we show the critical and caustic curves of 
the lens model overlaid the deconvolved R image. 



7. Discussion 

Extinction curves have been derived for the Galaxy and 
the Magellanic Clouds, mainly through multi-color pho- 
tometry of individual stars, for which the intrinsic lumi- 
nosities is known from their spectral types (Cardelli et al. 
19891). The extinction curves of the LMC and SMC can 



be fitted by the Galactic extinction law (see Mathis 1990; 
Fitzpatrick 1999). The main difference is the 2175 A bump 
(cf. Fig. ||) which is not observed in the LMC and SMC. 
At distances > 10 Mpc study of individual stars is diffi- 



cult and other methods must be considered. In this paper 
we have shown that accurate extinction curves of high- 
z galaxies can be determined by studying gravitationally 
lensed systems. We have found that the extinction curve 
of B1152+199 is well fitted by the Galactic extinction law. 
It is interesting to note that the inferred 1.3 < Ry < 2.1 
in principle is outside the validity of the parameteriza- 
tion of Galactic extinction curve by Cardelli et al. (1989). 



Note also that the differential extinction of Bl 152+199 
0.9 < E(B — V) < 1.1 , is one of the largest observed to 



date (Falco et al. |l999|) . 

The evidence for a Galactic extinction curve is intrigu- 
ing, but U and B band observations with longer exposure 
times and better seeing are needed to rule out other ex- 
tinction curves. These will sample the restframe UV ra- 
diation which is highly sensitive to dust extinction, and 
thus strongly constrain the shape of the extinction curve, 
and allow for the possible detection of the 2175 A bump. 

The predicted time delay of 59 days makes B1152+199 
an interesting candidate for monitoring. 
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